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This  report  was  prepared  as  part  of  USDI  Bureau  of  Land  Management  RDTS  No. 
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INTRODUCTION 


Streamflow  response  to  rainfall    is  controlled  by  the  net  water  budget  and  the 
routing  of  rainfall   to  stream  channels.     In  rain-dominated  forested  regions, 
including  the  Oregon  Coast  Range,  water  budget  aspects  of  forest  hydrology 
are  better  quantified  than  water  routing.     As  a  result,  most  of  the  commonly 
used  forest  hydrology  models  have  been  designed  to  assist  in  the  analysis  of 
seasonal   or  annual   runoff  volumes  (Goldstein,  et  al .  1974;   Troendle  and  Leaf, 
1980;   Silvey  and  Rosgen,  1981).     The  effects  of  timber  harvesting  on  runoff 
volumes  are  reflected  primarily  through  reductions  in  evapotranspi rational 
demand  and  the  corresponding  increases  in  soil  water  content. 

An  important  management  concern  is  the  effect  of  timber  harvesting  on  the 
size  of  large  peak  flows  and  the  duration  of  flows  greater  than  the  bankfull 
discharge.     There  is  a  growing  body  of  literature  which  suggests  that  it  is 
not  simply  runoff  volume,  but  flows  which  exceed  the  bankfull   channel 
capacity  which  most  strongly  influence  bed  material   transport,  channel 
morphology  and  channel   stability  (Emmett,  1976;   Dunne  and  Leopold,  1978; 
Jackson  and  Beschta,  1982;   Andrews,  1983).     Other  discharge-related 
management  issues  for  small   forest  streams  include  prediction  of  design  flows 
for  culvert  and  bridge  design,   floodplain  delineation,  and  low  flow  analysis 
for  water  supply  and  fisheries. 

There  are  no  widely  used,  physically  based  storm-period  runoff-routing  models 
for  the  Oregon  Coast  Range.     Thus,  certain  consequences  of  timber 
harvesting,   such  as  roads  and  skid  trails,  which  change  the  infiltration 
properties  of  forest  soils  and  which  may  influence  the  routing  of  rainfall   to 
streams,  are  difficult  to  analyze  in  terms  of  their  hydrologic  effects. 

Furthermore,  the  effects  of  timber  harvest  on  other  storm-period  hydrologic 
variables  such  as  individual   peak  flows  and  the  duration  of  very  high  flows 
are  not  readily  predicted.     Altered  hydrologic  regimes  in  turn  affect  channel 
processes  and  influence  instream  biologic  values. 

The  purposes  of  the  report  are  to: 

1)         Describe  the  rainfall-runoff  processes  in  the  Oregon  Coast  Range, 
and  review  the  results  of  studies  in  Western  Oregon  on  the  effects 
of  timber  harvesting  on  runoff. 

Review  and  analyze  available  techniques  for  predicting  the 
hydrologic  effects  of  timber  harvesting  in  rainfall -dominated 
regions. 


2) 


3)  Provide  a  review  of  the  issue  of  sediment  routing,   channel 
stability,  and  the  relationship  between  hydrology  and  channel 
morphology. 

4)  Provide  additional   analysis  of  the  Alsea  Watershed  Study  data 
(Harris,  1977)  to  better  describe  the  effects  of  timber  harvesting 
on  high-flow  durations  and  assess  the  potential   for  developing  and 
validating  simple  storm-period  rainfall -runoff  prediction  procedures, 


We  conclude  that  available  techniques  adequately  predict  the  effect  of  tim- 
ber harvest  on  annual  and  seasonal  water  yields.  Techniques  for  estimating 
design  flows  for  instream  structures  (e.g.  culverts,  bridges)  are  less  well 
developed,  but  generally  adequate  for  most  forest  hydraulic  engineering 
purposes.  Techniques  for  predicting  the  effects  of  timber  harvest  on  peak 
flows  or  channel -modifying  flows  (defined  as  flows  greater  than  the  bankfull 
discharge)  are  not  well  developed.  Furthermore,  there  are  no  subjective  or 
deterministic  procedures  for  quantifying  changes  in  channel  form  or  function 
given  changes  in  channel -forming  flows.  It  may  be  difficult  to  separate  the 
effects  of  altered  hydro! ogic  regime  from  changes  in  sediment  delivery  and 
sediment  transport  when  analyzing  channel  responses  to  timber  harvest 
activities. 


BACKGROUND  LITERATURE  REVIEW 


Coast  Range  Hydrology 


The  hydrology  of  small,  low-elevation  forest  streams  in  the  Oregon  Coast 
Range  is  dominated  by  low- to-moderate  intensity,  long-duration  frontal 
rainstorms,  and  shallow  subsurface  transmission  of  water  to  stream  chan- 
nels. A  descriptive  model  of  the  subsurface  stormflow  process  was  developed 
by  Hewlett  and  Hibbert  (1967)  and  further  described  by  Dunne  (1978)  and  Abdul 
and  Gillham  (1984).  In  that  model,  streamflow  response  to  precipitation 
occurs  rapidly  through  a  process  called  transelatory  flow.  As  rainfall  over 
a  watershed  increases,  the  area  directly  contributing  to  streamflow,  the 
source  area,  increases  and  the  active  channel  network  expands  in  length  and 
width.  The  expansion  and  contraction  of  drainage  networks,  and  stream  dis- 
charge rates  are  controlled  by  precipitation  characteristics,  the  amount  of 
water  stored  in  forest  soils,  and  subsurface  flow  characteristics  of  the 
watershed.  In  general,  as  the  amount  of  water  stored  in  forest  soils 
increases,  the  greater  is  the  percentage  of  rainfall  for  any  given  storm 
which  appears  as  stormflow  runoff  (quickflow),  Streamflow  response  to 
precipitation  is  directly  related  to  slope  steepness  and  drainage  density 
(Harr,  1983(a)). 

About  75  percent  of  the  annual  precipitation  and  runoff  in  the  Oregon  Coast 
Range  occurs  between  October  1  and  April  1.  Average  annual  precipitation  is 
strongly  influenced  by  elevation  and  ranges  from  40  in.  on  the  east  edge  of 
the  Coast  Range  to  over  100  in.  at  higher  elevations.  Snowfall  does  occur, 
but  is  generally  transitory  at  lower  elevations.  Melting  snowpack  during 
rainfall  has  contributed  to  some  of  the  largest  peak  flows  in  the  region. 
Annual  runoff  is  roughly  70  percent  of  annual  precipitation.  Most  rainfall 
which  does  not  appear  as  runoff  is  lost  to  interceDtion  (roughly  8-20 
in./yr.)  and  evapbtranspiration  (roughly  6-20  in. /yr.)  (Harr,  1983(a)).  Both 
interception  and  evapotranspiration  losses  increase  with  increased  canopy 
cover  and  increased  annual  precipitation. 

Harris,  et  al .  (1979)  classified  the  entire  Oregon  Coast  Range  as  a  single 
hydrologic  region  for  purposes  of  flood-frequency  analysis  (Figure  1).  The 
region  has  similar  climatic  and  topographic  characteristics.  "Geologic  parent 
materials  consist  primarily  of  sedimentary  rocks  of  the  Tyee  and,  to  a  lesser 
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extent,  the  Umpqua,  Yamhill   and  Nestucca  formations.     There  are  also 
scattered  areas  of  volcanic  igneous  rocks  and,   in  the  extreme  southern 
portion  of  the  range,  considerable  metamorphic  geologic  material.     Soils  are 
generally  deep  and  range  in  texture  from  clay  loams  to  sandy  and  gravelly 
loams.     The  Coast  Range  is  sharply  dissected  and  has  high  drainage  densi- 
ties, which  result  in  rapid  movement  of  water  to  stream  channels. 

First  order  stream  channels  are  usually  steep  and  incised  to  bedrock.     Second 
and  third  order  streams  are  often  controlled  structurally  by  large  logs, 
which  serve  to  trap  fluvial   sediments  and  dissipate  stream  energy  in  steps 
(Swanson  and  Lienkaemper,  1978).     These  channels  are  often  adjustable  and 
dynamic  during  periods  of  high  runoff. 

Campbell   and  Sidel   (1984)  analyzed  steam  gage  records  in  the  coast  hydro- 
logic  region  of  Oregon  (Fig.   1)  as  well   as  5  other  hydrologic  regions  and 
found  the  log-Pearson  type  III  distribution  with  regional   skew  coefficients 
to  be  suitable  for  use  in  all   regions  tested.     The  following  flood-frequency 
regression  equations  were  developed  for  the  Coast  Region  (Campbell  and  Sidel, 
1984). 

Q10  =  O.lllA1-0^'-49 
Q25  =  0  ^a"*  -01  £.51 
Q50  =  o! 152A1-01p.50 
QlOO  =  0J66A»'00E-50 


where 

Q  =  flow,  m2/s 

A  =  drainage  area,  km2 

E  =  elevation,  m 

Lystrom  (1970)  also  developed  multiple  regression  equations  for  different 
flood  peaks  for  gaged  rivers  in  Oregon  based  upon  watershed  characteristics 
(e.g.,  slope,  length,  storage,  elevation,  area,  precipitation,  precipitation 
intensity,  etc. ). 

Complete  reviews  of  the  hydrology  of  western  Oregon  forested  regions  are 
provided  in  Harr  (1976(b))  and  Harr  (1983(a)). 

Effects  of  Logging  Activities  on  Hydrology 

Logging  activities  -  timber  harvesting,  road  building  and  skidding  -  may 
affect  almost  all  aspects  of  the  forest  hydrologic  cycle  including  net 
precipitation,  infiltration,  evaporation,  transpiration  and  the  delivery  of 
rainfall  to  streams.  Two  research  programs,  the  Al sea  Watershed  Studies  and 
the  H.  J.  Andrews  Experimental  Forest,  have  provided  most  of  the  quantified 
information  on  the  hydrologic  response  to  logging  and  road  building  in  the 
Douglas  Fir  zone  of  western  Oregon.  The  Alsea  watersheds,  located  in  the 
Coast  Range  near  Newport,  Oregon,  provide  data  most  representative  of  the 
Coastal  hydrologic  region. 
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A  brief  summary  is  provided  below  on  the  effects  of  logging  on  specific 
aspects  of  the  hydrologic  cycle:  interception,  evapotranspiration,  and 
infiltration   Streamflow  response  to  these  changes,  including  annual  water 
yield,  peak  flows  and  low  flows  are  then  reviewed.  The  discussion  is  limi- 
ted to  the  rainfall -runoff  situation  because  of  its  predominance  in  low 
elevation  coastal  watersheds. 

Interception 

The  amount  of  water  intercepted  by  vegetation  cover  is  a  function  of  the 
form,  density  and  surface  texture  of  the  vegetation.  The  amount  of  precipi- 
tation lost  through  interception  is,  in  turn,  equal  to  the  amount  of  water 
evaporated  from  interception  storage.  .Storm  period  interception  loss  in- 
creases with  rainfall  duration,  but  the  proportion  of  rainfall  lost  to 
interception  decreases  as  rainfall  amounts  increase  over  the  storage  capactiy 
of  the  vegetation.  Annual  interception  loss  increases  with  rainfall 
frequency  and  losses  are  greatest  at  the  beginning  of  storms. 

A  common  expression  for  storm  period  interception  loss,  I,  is 

I  =  c(l-e-P/c)  +  VEt  (2) 

where: 

c  =  interception  storage  capacity 

p  =  precipitation  amount 

V  =  number  of  vegetation  surfaces 

E  =  evaporation  rate 

t  =  time 

In  the  Coast  Range,  c  varies  from  less  than  .01  inches  for  deciduous  trees 
without  foliage  to  roughly  .05  inches  for  mature  Douglas  fir.  Clearcut 
logging  reduces  c  to  essentially  zero.  The  recovery  of  c  to  mature  forest 
levels  is  related  to  species  composition  and  the  rate  of  the  reforestation 
process. 

Increases  in  net  precipitation  following  logging  due  to  interception  may  be 
partially  offset  by  reduced  fog  drip  (Harr,  1982). 

Evapotranspi  rati  on 

Evapotranspiration  is  the  process  by  which  plants  consume  water  from  the  soil 
and,  ultimately,  vaporize  it  to  the  atmosphere.  The  potential  evapo- 
transpirational  rate,  Es,  is  a  maximum  loss  rate  unlimited  by  soil  or  plant 
factors,  but  controlled  by  meterologic/energy  factors.  Actual  trans- 
piration, Ea,  becomes  less  than  Es  when  water  availability  in  the  soil 
becomes  limiting.  Leaf  and  Brink  (1975)  consider  that  Ea  =  Es  when  avail- 
able soil  water  is  more  than  50  percent  of  field  capability.  The  ratio  Ea/Es 
then  decreases  linearly  to  zero  as  available  soil  water  is  reduced  to  the 
wilting  point. 

In  the  Alsea  watersheds  of  the  Coast  Range,  Harr  (1983)  estimated  actual 
annual  evapotranspiration  was  12  to  18  inches  per  year.  In  oTd  "growth 
Douglas  fir  forest  at  the  H.  J.  Andrews  forest,  Rothacher  (1963)  estimated 


annual  evapotranspiration  was  roughly  16  to  18  Inches  per  year.  Clearcut 
logging  reduces  Ea  to  approximately  zero  the  first  year  following  logging 
(some  bare  soil  evaporation  will  occur).  The  results  of  both  the  Alsea  study 
and  watershed  studies  at  H.  J.  Andrews  and  elsewhere  suggest  that  Ea  recovers 
rapidly  to  pre-logging  levels  following  reforestation.  While  recovery  is  a 
function  of  species  type,  site  conditions  and  subsequent  management,  Ea  may 
approach  pre-logging  levels  in  roughly  5  to  15  years  in  western  Oregon 
Douglas  fir  forests. 

Infiltration 

Infiltration  capacities  in  undisturbed  Coast  Range  forests  are  high. 

Saturated  hydraulic  conductivities  in  excess  of  118  inches  per  hour  have  been 

measured  in  Oregon  Coast  Range  forests  (Yee,  1975),  explaining  the  general 
absence  of  overland  runoff. 

Road  building  and  tractor  skidding  cause  compacted  soils,  a  loss  of  surface 
soil  structure  and  reduction  of  macro  pores.  This  results  in  greatly  re- 
duced infiltration  capacities  and  increases  the  possibility  for  overland 
runoff  (Johnson  and  Beschta,  1980).  In  some  soils,  the  burning  of  slash 
following  logging  may  also  reduce  infiltration  capacities  by  creating  a 
temporary  hydrophobic  condition.  When  soil  disturbance  results  in  surface 
runoff,  surface  erosion  occurs  and  the  translation  of  rainfall  to  streamflow 
may  occur  more  quickly  -  possibly  affecting  instantaneous  flow  rates. 

Surface  runoff  rates  are  usually  high  from  road  surfaces  (Reid  and  Dunne, 
1984).  However,  the  runoff  is  susceptible  to  management  through  road 
placement,  design,  and  drainage  control;  thus  the  effects  on  streamflow  are 
^ery   difficult  to  generalize. 

The  reduction  in  infiltration  on  skid  trails  is  a  function  of  the  type  of 
yarding  (e.  g.,  tractor,  cable)  -  which  affects  the  amount  of  compacting 
energy  translated  to  the  soil  -  soil  texture,  water  content  and  organic 
matter  content.  Generally,  tractor  skidding  on  moist,  fine-textured  soils 
causes  the  greatest  reduction  in  infiltration.  Recovery  to  pre-logging 
infiltration  capacities  is  dependent  upon  revegetation  success  and  other 
soil-influencing  processes  such  as  freeze-thaw.  Rehabilitation  of  infil- 
tration capacities  on  skid  trails  is  also  subject  to  management  through 
mechanical  tillage  or  ripping,  so  effects  on  stream  flow  can  be  managed. 

Annual  Water  Yield 

The  effects  of  logging  on  annual  and  seasonal  water  yield  in  the  Douglas  fir 
region  of  western  Oregon  have  been  well  summarized  by  Rothacher  (1971),  Harr 
(1976(a)),  Harr  (1983(a))  and  Harr  (1983b),  and  are  only  briefly  reviewed 
here. 

Annual  water  yield  increases  at  four  clearcut  watersheds  at  H.  J.  Andrews  and 
Alsea  watersheds  ranged  from  about  4  inches  per  year  to  over  23  inches  per 
year.  Eighty-seven  percent  of  variation  in  the  annual  increases,  Y,  at  the 
H.  J.  Andrews  studies  was  explained  with  a  two  parameter  linear  model 


incorporating  number  of  years  after  logging,   X] ,  and  annual   precipitation, 
x2- 

Y  =  31.41    -  2.08X]    +  0.091    X2  (3) 

In  general,  water  yield  increases  are  greatest  in  wetter  years  and  in  the 
fall   and  spring,  and  decrease  over  time  following  the  recovery  of  vegeta- 
tion.    Reductions  in  interception  and  evapotranspiration  losses,  as  de- 
scribed above,  account  for  the  observed  increases  in  yield. 

We  developed  a  similar  regression  model   for  annual   flow  increases  at  Needle 
Branch  in  the  Alsea  Watersheds  following  logging: 

Y  =     -0.32  -  0.21X]    +  0.17  X2  (4) 

r2  =  0.42 
(Std.   error  of  the  estimate  of  Y  =  3.8  in.) 

However,  the  post-logging  period  was  only  7  years  at  Needle  Branch  which  is 
probably  inadequate  for  this  type  of  analysis. 

Increases  in  annual  water  yields  have  also  been  measured  for  partially  logged 
watersheds  (Harr,  1983b).  Increases  occurred  in  both  patch-cut  and 
shelterwood-cut  watersheds,  but  were  generally  less  than  for  clearcut  water- 
sheds. While  increases  in  water  yield  can  be  expected  to  be  related  to  the 
percent  of  the  watershed  logged,  significant  increases  in  annual  water 
yields  were  not  detected  at  one  watershed  at  the  H.  J.  Andrews  forest  until 
40  percent  of  the  watershed  was  clearcut.  When  a  smaller  percent  of  the 
watershed  is  logged,  it  is  possible  that  reductions  in  Es  are,  in  part, 
compensated  for  by  increases  in  Ea  in  the  remaining  stand.  Also,  the  loca- 
tion of  a  cut  within  a  watershed  may  also  influence  its  effect  on  runoff. 

Peak  Flows 

The  effects  of  logging  on  the  size  of  peak  flows  -  particularly  large  peak 
flows  -  are  not  easily  generalized.  Changes  in  seasonal  or  annual  runoff 
volumes  following  timber  harvest  are  well -documented  and  well -explained  by 
net  changes  in  the  water  budget.  However,  the  timing,  or  rate,  of  runoff  is 
another  variable  which  influences  the  shape  of  the  hydrograph  and  the  size  of 
peak  flows.  Therefore,  to  predict  the  effects  of  logging  on  peak  (or 
instantaneous)  flows,  it  is  necessary  to  : 

(1)  Describe  the  effects  of  an  altered  water  budget  on  the  volume  of 
storm-period  runoff,  and 

(2)  Describe  the  effects  on  the  routing  of  rainfall  to  stream 
channels  (i.e.,  the  distribution  of  runoff  over  time). 

The  results  of  existing  watershed  studies  strongly  suggest  that  the  effects 
of  logging  on  peak  flows  in  this  region  are  to: 


(1)  Increase  small  peak  flows  in  early  fall  and  spring  due  to  increased 
runoff  volumes  which  are  the  result  of  higher  soil  water  contents  in 
logged  areas  than  forested  areas  during  those  seasons,  and 

(2)  Increase  larger  peak  flows  which  result  from  rainfall  on  saturated 
soils  in  those  situations  where  the  routing  of  rainfall  to  streams 
has  been  accelerated  (Harr,  et  al .  1975;  Harr,  1976(a)). 

In  addition,  there  is  a  growing  body  of  literature  which  suggests  that  the 
snow  accummulation  and  melt  process  (including  rainfall  on  snow)  is  suffi- 
ciently affected  by  logging  to  influence  the  size  of  peak  flows  in  this 
region  (Harr  et  al .  1982,  Christner  and  Harr,  1982;  Harr  and  Berris,  1983). 

The  principal  mechanism  for  increasing  the  rate  of  delivery  of  rainfall  to 
stream  channels  is  to  reduce  infiltration  to  the  point  where  overland  runoff 
occurs  -  thus  altering,  in  part,  the  routing  mechanism  from  interflow  to 
overland  flow.  This  occurs  primarily  on  road  surfaces  and  severely  com- 
pacted soils  such  as  may  occur  on  skid  trails  and  landings. 

In  those  cases  where  large  peak  flows  have  been  significantly  increased 
(e.g.,  A! sea  study  -  Deer  Creek  3),  large  proportions  of  the  watershed 
(  12%)  were  occupied  by  roads,  skidtrails  and  landings  (Harr,  1975;  Harr, 
1976(a);  Harr,  1983(a)).  In  studies  of  the  effects  of  cutting  alone,  when 
roads,  trails  and  landings  occupied  a  negligible  percent  of  the  watershed, 
large  peak  flows  were  generally  unaffected,  even  though  smaller  peak  flows 
(and  the  average  of  all  peak  flows)  were  increased  (Rothacher,  1971;  Harris, 
1977;  Harr  et  al . ,  1975;  Harr,  1976(a);  Harr,  1983(a)). 

Whether  roads,  skidtrails  and  landings  increase  large  peak  flows  and  how 
large  any  increases  will  be  is  difficult  to  predict  and  is  dependent  upon 
quantification  of  both  the  pre-  and  post-logging  hydro! ogic  routing  systems 
and  percent  of  total  watershed  area  these  facilities  occupy.  While  models  of 
both  subsurface  storm  flow  routing  (Sloan,  et  al .  1983)  and  surface  runoff 
(e.  g.,  Simons,  et  al .  1975)  exist,  they  have  not  been  combined  and  validated 
in  this  region.  It  is  probable,  however,  that  the  percent  area  in  roads, 
skid  trails  and  landings  is  only  one  variable  to  be  considered  and  that  the 
location,  design  (including  drainage  design)  and,  if  appropriate,  reclamation 
of  these  features  are  also  important  factors  influencing  routing.  All  of 
these  factors  are  subject  to  management.  Thus,  effects  of  timber  harvest  on 
peak  flows  may  be  very   difficult  to  generalize  or  quantify. 

Low  Flows 

The  effects  of  logging  on  low  flows  are  summarized  in  Harr,   (1976(a)),   Harr 
(1983(a))  and  Harr  (1983b).     Increased  low  flows  (or  decreases  in  the  number 
of  low  flow  days)  occurred  following  logging  at  the  H.  J.  Andrews  and  the 
Alsea  watersheds.     Similar  increases  have  been  measured  at  other  sites 
in  western  Oregon.     While  the  increases  in  low  flows  were  proportionally 
large  (roughly  100%)  between  treated  and  controlled  watersheds,  they  ac- 
counted for  only  a  small   amount  of  the  total   annual  water  yieTd  increases. 


Furthermore,  the  increase  were  short  lived,  with  recovery  to  pre-logging 
conditions  occurring  in  roughly  five  years. 

Increases  in  low  flows  are  generally  explained  by  reduced  evapotranspi ra- 
tional demand  following  logging.  However,  in  at  least  one  case  at  the  Bull 
Run  principal  watershed  near  Portland,  Oregon,  low  flows  were  decreased 
following  patch  cutting.  This  decrease  has  been  attributed  to  reduced  pre- 
cipitation in  the  form  of  fog  drip  (Harr,  1980;  Harr,  1982).  Revegation  of 
streamsides  by  phreatophytes  may  also  reduce  low  flows  compared  to  the  pre- 
logged  condition  until  after  they  become  overtopped  by  Douglas  fir  and 
Hemlock  (Harr,  1983(b)). 

Effects  of  Hydrology  on  Channel  Processes 

Flow,  Qw,  and  sediment  discharge,  Qs,  are  the  primary  variables  influencing 
the  morphology  of  stream  channels  with  adjustable  beds  and  banks  (Schumm, 
1971).  Alluvial  channels  are  considered  to  be  "graded"  if  their  resulting 
hydraulic  characteristics  are  the  minimum  required  to  transport  the  sediment 
delivered  to  the  channel  (Leopold  and  Maddock,  1953).  Channel  geometry  and 
hydraulic  characteristics  will  change  when  either  flow,  sediment  discharge  or 
both  increase.  Schumm  (1971)  proposes  that  a  channel  will  widen  when  either 
Qw,  Qs  or  both  increase.  Channel  depth  will  increase  with  an  increase  in  Qw 
but  will  decrease  when  Qs  increases. 

Based  upon  the  above  criteria,  a  stream  channel  should  respond  to  an  in- 
creased flow  regime  by  becoming  wider  and  deeper.  However,  increased  flows 
following  logging  are  almost  always  accompanied  by  increased  sediment 
delivery  to  streams  (Brown,  1980).  Thus,  in  most  documented  changes  in 
channel  characteristics  following  logging,  the  effects  of  increased  Qs 
predominate  and  channels  tend  to  become  wider  and  shallower  (Lyons  and 
Beschta,  1983;  Beschta,  1984).  Jackson  and  Beschta  (1984)  suggest  that  the 
initial  response  of  a  stream  channel  to  increases  in  Qs  will  be  to  increase 
its  hydraulic  efficiency  by  reducing  form  roughness.  Reduced  form  roughness 
is  accomplished  by  reducing,  or  smoothing,  the  channel's  pool-riffle  pat- 
tern. This  permits  increased  velocities  and  an  increased  capacity  to 
transport  water  and  sediment.  Subsequent  responses  are  reduced  channel 
stability,  increased  channel  width,  and  decreased  channel  depth.  Reduction 
or  elimination  of  pool -riffle  patterns,  reduced  bed  and  bank  stability  and 
changes  in  particle  size  characteristics  of  substrate  all  influence  the 
biologic  quality  of  the  channel  (Moring  and  Lantz,  1975). 

While  the  above-channel  responses  may  occur  following  logging  if  there  are 
sufficient  increases  in  Qw  and  Qs,  the  hydraulic  geometry  equations  cannot  be 
solved  explicitly  -  therefore  the  exact  response  is  difficult  to  predict.  It 
becomes  even  more  difficult  to  separate  out  the  effects  of  increased  flows, 
alone.  We  do  know,  however,  that  not  all  flows  affect  channel  processes,  and 
therefore,  channel  morphology.  There  are  flows  below  which  there  is  \/ery 
little,  if  any,  movement  of  bed  or  bank  material.  When  a  channel  has  a  " 
sufficient  supply  of  large  bed  material  armoring  occurs  and  flows  required  to 
initiate  bed  material  movement  (critical  flows)  may  become  quite  "large." 
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Several  studies  have  provided  data  which  suggest  that  critical  flows  for  many 
alluvial  channels  approximate  the  bankfull  discharge  and  have  roughly  a  1.5 
to  2-year  return  period.  Andrews  (1983),  in  a  study  of  24  self- formed 
gravel -bed  rivers  in  Colorado,  determined  the  critical  dimension! ess  shear 
stress  required  to  entrain  given  particle  sizes.  He  found  that  while  smal- 
ler particles  were  frequently  entrained,  the  larger  (e.g.,  dgg)  particle 
sizes  were  not  entrained  until  roughly  the  bankfull  discharge.  At  the  Snake 
and  Clearwater  Rivers  in  Idaho,  Emmett  (1976)  observed  that  bed-load  trans- 
port rate  was  a  different  (and  higher)  function  of  unit  stream  power  when  the 
stream  was  capable  of  transporting  all  particle  sizes  than  at  lower  stream 
powers  where  bed  armoring  progressively  limited  the  bed  material  available 
for  transport.  Sands  were  transported  at  lower  flows,  whereas  coarse  gravels 
were  primarily  transported  upon  disruption  of  the  armor. 

Jackson  and  Beschta  (1982)  described  bed  material  routing  in  a  sand  and 
gravel -bed  channel  at  Flynn  Creek  in  the  Alsea  watersheds.  At  discharges 
less  than  bankfull,  bedload  transport  consisted  primarily  of  sand-sized 
materials  being  routed  over  stable  gravel  riffles.  At  approximately  the 
bankfull  discharge,  riffle  materials  were  also  transported  downstream. 
However,  because  of  nonuniform  channel  geometries,  bedload  transport  rates 
were  unsteady  and  several  sequences  of  partial  riffle  scour  and  fill  oc- 
curred. Upon  the  recession  of  high  flows,  a  final  riffle  deposition  and  a 
natural  sorting  of  particle  sizes  occurred.  The  bed  material  composition  of 
the  riffles  was  considerably  more  coarse  than  that  of  the  bedload  sediment  in 
transport  at  the  time  riffle  sediments  are  deposited.  "Left-over"  sand- 
sized  bed  material  therefore  deposited  elsewhere  in  the  stream  channel  at 
lower  flows  -  generally  in  pools,  channel  edges  and  backwaters. 

In  summary,  increases  in  Qs  and  Qw  seem  to  be  the  principal  factors  effect- 
ing changes  in  channel  morphology.  While  these  two  factors  interact  to  in- 
fluence channel  geometry  and  channel  processes,  net  channel  response  seems  to 
reflect  the  effect  of  increased  Qs  over  increased  Qw.  There  is  evidence  that 
channels  are  most  responsive  to  flows  greater  than  approximately  bank-  full, 
when  stream  powers  are  sufficient  to  mobilize  the  entire  range  of  sizes  of 
bed  and  bank  materials.  Lower  flows  serve  primarily  to  transport  sand  sized 
materials  and  redistribute  sands  to  pools  and  other  areas  of  relatively  low 
shear  stresses. 


11 


REVIEW  AND  ANALYSIS  OF  QUANTITATIVE  TECHNIQUES 
FOR  ASSESSING  HYDROLOGIC  EFFECTS  OF  LOGGING  ACTIVITIES 
Introduction 

A  comprehensive  procedure  for  evaluating  the  effects  of  si Ivi cultural  activ- 
ities on  the  hydrologic  cycle  is  provided  in  WRENS  (Troendle  and  Leaf, 
1980).  WRENS  may  be  used  to  analyze  changes  in  the  seasonal  or  annual  water 
budget  for  forested  watersheds  and  can  produce,  as  output,  change  in  the 
6-day  flow  duration  curve  for  a  watershed.  The  hydrology  component  of  WRENS 
uses  two  existing  models.  The  PROSPER  model  (Goldstein,  et  al . ,  1974)  is 
used  for  rainfall -runoff  situations  and  the  WATBAL  model  (Leaf  and  Brink, 
1973)  is  used  for  snowmel t-runoff  situations.  Other  procedures,  for  example 
the  HYSED,  procedure  developed  for  Region  2  of  the  Forest  Service  (Silvey  and 
Rosgen,  1981)  use  modified  versions  of  WRENS  for  their  hydrologic  com- 
ponents. Most  procedures  incorporating  routing  routines  which  are  capable  of 
analyzing  individual  storm  events  and  instantaneous  flows  are  either  still 
under  development,  or  are  not  validated  for  western  Oregon.  Sloan  et  al . 
(1983)  evaluated  several  subsurface  -  storm-period  models  on  steep-  sloped 
forested  watersheds  in  the  eastern  United  States.  While  their  models  have 
not  been  validated  in  Western  Oregon,  they  may  provide  a  basis  for  eventual 
development  of  a  physically  based,  storm-period  model  for  this  region. 
Overton  and  White  (1977)  developed  a  storm-period  runoff  model  for  two 
watersheds  in  the  H.  J.  Andrews  forest,  but  discontinued  model  refine-  ment 
due  to  a  lack  of  information  on  several  ecosystem  processes  and  a  need  for  a 
snowmel t  subroutine. 

One  method  was  reviewed  which  predicted  changes  in  water  yield  and  channel 
stability  for  the  Kootenai  National  Forest  (Galbraith,  1973).  This  method, 
however,  assumes  that  water  yield  -  not  necessarily  flows  over  some  level  - 
affect  channel  stability.  This  assumption  is  not  well -supported  in  the 
literature. A  more  detailed  review  of  the  WRENS  procedure,  in  particular  the 
PROSPER  model,  is  provided  below.  Additional  reviews  of  the  Overton  and 
White  model  and  the  models  reported  in  Sloan,  et  al .  are  also  provided. 

WRENS 

The  hydrology  component  of  WRENS  was  developed  to  assist  professionals  in  the 
analysis  of  the  effects  of  logging  -  primarily  vegetation  removal   -  on 
seasonal   and  annual   runoff  volumes.     Model   parameters  have  been  generalized 
for  eight  physiographic  regions  -  one  of  which  encompasses  the  coastal 
forests  of  northern  California,  Oregon  and  Washington.     For  rain-dominated 
regions  below  3,000  -  4,000  feet  in  the  Pacific  Coast,  the  PROSPER  model 
(Goldstein,  et  al . ,  1974)  is  applied.     The  model   has  been  tested  on  Needle 
Branch  of  the  Alsea  Watersheds,  with  the  objective  of  extrapolating  local 
observations  for  regional   use. 

PROSPER  is  strictly  a  water  balance  model  with  no  surface  or  subsurface 
routing  components.     PROSPER  emphasizes  the  evapotranspiration_  component  of 
the  hydrologic  cycle  by  integrating  the  interaction  between  available  soil 
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water,  precipitation,  and  energy.     Regional   evapotranspiration  data  is  used 
by  the  model.     Leaf  area  index  is  used  as  the  primary  evapotranspiration 
modifier.     In  addition,  the  evapotranspiration  function  is  affected  by 
changes  in  available  water  caused  by  altering  rooting  depth.     Water  avail- 
able for  streamflow  is  determined. by  subtracting  seasonal   evapotranspiration 
from  net  precipitation.     The  model   can  handle  several  si Ivi cultural   pre- 
scriptions within  a  watershed. 

A  schematic  of  the  PROSPER  model   as  provided  in  the  WRENS  documentation 
(Troendle  and  Leaf,  1980)  is  reproduced  in  Appendix  I. 

The  WRENS-PROSPER  procedure  provides  a  simple  physically  based  procedure  for 
calculating  the  effects  of  logging  on  seasonal   and  annual   runoff.     It  has 
been  validated  at  Needle  Branch  in  the  Oregon  Coast  Range.     Because  regional 
evapotranspiration  potential   is  fairly  uniform  and  not  particularly  sensi- 
tive to  aspect  and  slope  during  the  rain  season,  and  because  it  can  handle 
differences  in  watershed  soils,  WRENS-PROSPER  should  be  fairly  transferable 
to  other  watersheds  in  the  Coast  Range.     Functions  relating  evapotranspira- 
tion to  leaf  area  index  were  developed  in  the  southeastern  United  States. 
They  still   need  additional   validation  in  the  Oregon  Coast  Range. 

The  main  drawback  to  WRENS-PROSPER  is  that  the  model   output  is  not  particu- 
larly useful   in  the  analysis  of  many  hydrology  related  management  issues 
including  design-flow  hydrology,  floodplain  hydrology,  or  channel   processes. 

Overton-White 

Overton  and  White  (1978)  developed  several   versions  of  a  physically-based 
storm-period  rainfall -runoff  model    for  the  H.  J.  Andrews  experimental   forest 
as  part  of  a  study  of  modelling  strategy.     Their  objective  was  to  develop  a 
submodel   of  an  ecosystem  model   and  as  such  were  less  concerned  with  the  ac- 
curacy of  the  hydro! ogic  output  than  with  the  capability  of  the  model   to 
interact  realistically  with  other  ecosystem  components.     The  model  was 
developed  for  Watershed  10  of  H.  J.  Andrews  Forest,  but  some  input  data  were 
generalized  from  Watershed  2  and  other  local  data  sets.     As  in  PROSPER,  the 
model    is  a  water  balance-type  model,  but  unlike  PROSPER  it  does  attempt  to 
represent  both  soil  water  flow  and  groundwater  lateral   flow.     Also,  by 
incorporating  an  infiltration  routine,  the  model    is  capable  of  generating  an 
overland  flow  component.     By  compartmentalizing  the  model   into  above-ground 
and  below-ground  components,   it  was  shown  that  the  above  ground  component  - 
primarily  evapotranspiration-controlled  runoff  volume,  and  the  below  ground 
component  -  which  generally  related  to  storage  routing  -  controlled  hydro- 
graph  timing  and  shape.     Model   output  is  daily  discharge  in  m3. 

Additional   descriptions  of  the  models  from  Overton  and  White  (1978)  are 
provided  in  Appendix  II. 

Overton  and  White  discontinued  model   refinement  because  they  felt  certain 
processes  required  more  complete  specification  and  because  they  felt  a  need 
for  an  improved  snowmelt  routine.     No  additional   attempts  have" been  found  in 
the  literature  to  further  validate  or  refine  the  Overton-White  models.     Also, 
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their  ability  to  predict  changes  due  to  logging  activities  have  not  been 
tested.     Expanding  source  areas  are  not  well   incorporated  conceptually  by  the 
Overton-White  models. 

In  addition  to  an  improved  snowmelt  routine,  additional  work  in  modelling 
interflow  and  a  capability  to  subdivide  watersheds  into  "routing"  or 
"response"  units  might  improve  the  Overton-White  models. 

Some  of  these  shortcomings  have  been  addressed  in  other  large,  distributed 
parameter  models  -  for  example  the  USGS-PRMS  modelling  system  (Leavensly  et 
al I.  1983).     These  models  also  have  fairly  primitive  subsurface  routing  rou- 
tines,  fail   to  represent  the  concept  of  expanding  source  areas  and  would 
require  calibration  and  validation  in  the  Coast  Range. 

Most  physically  based  distributed  parameter  storm-period  models  are  data 
intensive  and  not  conducive  to  quick  "planning  level"  hydrologic  analyses. 
However,  as  research  tools  and  when  applied  to  special   situations,  they  offer 
potential   for  far  more  detailed  insight  into  hydrologic  systems  -  including 
instantaneous  flows  -  than  simple  water  budget  models. 

Sloan,  Moore,  Coltharp  and  Eizel 

Sloan,  et  al .   (1983)   provide  a  thorough  report  on  modelling  surface  and 
subsurface  stormflow  on  steeply-sloping  forested  watersheds.     Their  ob- 
jective is  to  identify  and  develop  techniques  for  modelling  rainfall -runoff 
consistent  with  the  theories  of  variable  source  areas-interflow.     They  re- 
view the  Variable  Source  Area  Simulator  Model   (VSAS)  (Troendle  and  Hewlett, 
1979).     They  also  review  five  subsurface  flow  models  for  homogeneous  soils: 
two  of  the  models  were  based  on  the  Richards  equation,  on  on  the  kinematic 
wave  equation  and  two  storage  models  -  the  Boussenesq  and  kinematic  storage 
models.     Finally  the  authors  develop  and  test  a  simple  continuous  daily 
rainfall -runoff  model  based  on  the  theories  of  interflow  and  variable  source 
areas. 

VSAS 

The  Variable  Source  Area  Simulator  Model  (VSAS)  was  developed  for  small, 
forested  watersheds  and  is  based  on  the  concept  that  instantaneous  streamflow 
is  the  sum  of  subsurface  flow,  precipitation  on  channel  and  saturated  zones 
and  overland  flow  from  impervious  areas  (Troendle  and  Hewlett,  1979). 
Subsurface  flow  is  described  by  the  Darcy  and  Richards  equations,  where 
unsaturated  conductivity  is  related  to  water  content  by 

K(0)  =  a  e  b  (5) 

where  a  and  b  are  constants  (Green  and  Corey,  1971).     The  watershed  is 
partitioned  into  segments  and  elements  (including  at  least  two  soil   layers) 
and  flow  goes  from  one  element  to  another  based  upon  the  above  laws  of  sub- 
surface flow.     A  VSAS  model   required  detailed  descriptions  of  watershed  soils 
and  as  such  is  data-intensive  and  probably  best  suited  to  research. 
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Sloan,  et  al . 

A  simple,  daily  rainfall -runoff  model  for  forested  watersheds  was  developed 
by  Sloan,  et  al.  (1983).  The  model  consists  of  three  water  stores  -  an 
interception  zone,  a  soil  zone,  and  a  groundwater  zone.  All  net  precipi- 
tation is  assumed  to  enter  the  soil  zone  -  except  on  saturated  areas  where  it 
becomes  runoff.  The  size  of  the  saturated  zone  increases  as  the  soil  zone 
becomes  wetter.  Subsurface  drainage  from  the  soil  zone  increases 
exponentially  as  the  water  content  increases.  Evapotranspirational  loss  is 
limited  either  by  potential  evapotranspiration  or  available  soil  water.  A 
model  schematic,  function  descriptions  and  parameter  descriptions  are 
provided  in  Appendix  I. 

The  model  was  calibrated  and  tested  on  the  Little  Mi  11  seat  Watershed  in 
eastern  Kentucky  using  a  split  record  technique.  Daily  streamflows  were  well 
predicted  by  the  model. 

Application  of  the  model  requires  partitioning  a  watershed  into  a  series  of 
interconnected  water  stores,  and  as  such  is  somewhat  data  intensive.  How- 
ever, because  it  is  a  conceptual  lumped-parameter  representation  of  the 
rainfall -runoff  process,  it  has  easier  data  requirements  than  more  dis- 
tributed parameter  models  such  as  Overton-White  or  PRNIS.  Its  greatest 
strength  may  be  in  its  conceptual  handling  of  variable  source  areas.  The 
model  has  not  been  tested  in  western  Oregon  but,  because  of  its  represen- 
tation of  expanding  source  area-interflow  processes,  it  can  be  expected  to 
apply  to  the  Coast  Range  forested  situation.  Also,  its  ability  to  model  the 
effects  of  logging  has  not  been  tested.  By  adding  an  infiltration  routine 
(presently  all  net  precipitation  is  infiltrated)  the  model  should  be  suited 
to  an  analysis  of  the  effects  of  roads  and  skidtrails  on  runoff  because  it 
allows  the  watershed  to  be  subdivided  into  water  "stores."  Like  Overton- 
White  and  PROSPER,  the  model  needs  a  snowmelt  routine  and  channel  routing 
routine. 
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EFFECTS  OF  LOGGING  ACTIVITIES  ON  FLOWS  GREATER  THAN  BANKFULL: 
A  REVIEW  OF  THE  ALSEA  WATERSHED  STUDY 

Introduction 

The  key  variable  with  respect  to  high  flow  influences  on  channel  processes 
is  the  total  force  exerted  on  the  channel  bed  and  banks  in  excess  of  the 
critical  force  required  for  incipient  bed  or  bank  instability.  If  channel 
instability  occurs  at  bankful  or  greater  discharges,  we  can  derive  either 
average  bed  shear  stress,  T;    ,  or  unit  stream  power,  VS,  from  discharge  and 
channel  geometry  information  and  integrate  for  total  force  or  stream  power 
expended  during  bankful!  or  greater  flows.  If  we  assume  the  hydrograph  peak 
can  be  approximated  by  a  triangle,  synthetic  peak  flow  (greater  than  bank- 
full  flow)  hydrographs  can  be  developed  given  estimates  of  peak  flow,  qp, 
and^the  duration  of  bankful!  or  greater  flow,  Djj.  We  can  then  solve 
for  To    or  VS  given  relationships  between  discharge  and  average  flow  velocity 
or  discharge  and  hydraulic  radius. 

We  hypothesized  that  qn  should  be  predictable  from  precipitation  for  a 
saturated  watershed.  To  test  our  hypothesis,  we  tried  to  predict  peak  flows 
from  precipitation  and  antecedent  streamflow  on  the  alien  watersheds.  We 
also  hypothesized  that  DD,  would  also  be  predictable.  We  attempted  to 
test  the  effects  of  timber  harvest  on  DD  by  reviewing  available  data  from 
the  alien  watersheds.  Our  analyses  are  exploratory  and  not  complete. 
However,  we  wanted  to  determine  if  additional,  more  thorough  analylses  of 
the  type  might  improve  our  ability  to  quickly  estimate  qp,  and  Db,  and 
eventually  the  effects  of  timber  harvesting  on  channel  modifying  flows, 
without  having  to  apply  more  cumbersome  instantaneous  flow  models. 

Review  of  A! sea  Results 

The  effects  of  logging  on  high  flows  at  two  of  the  Alsea  watersheds  have 
been  well  analyzed  by  Harris  (1977)  and  Harr,  et  al .  (1975). 

Harris  analyzed  all  peak  flows  greater  than  50  csm  (roughly  bankful! 
discharge)  on  Flynn  Creek  (the  control)  and  concurrent  peaks  on  Needle 
Branch  and  Deer  Creek.  The  mean  of  these  larger  peaks  was  increased  by  17.8 
csm  over  the  predicted  mean  of  91.2  csm  on  Needle  Branch  following  82 
percent  clearcutting  and  burning.  This  increase  was  not  considered 
significant  (p  =  0.95).  Five  percent  of  the  Needle  Branch  was  in  roads  or 
skid  trails.  The  mean  large  peak  flow  on  Deer  Creek  was  increased  by  1.8 
csm  over  the  predicted  mean  of  78.5  csm  following  25  percent  patch  cutting. 
Four  percent  of  Deer  Creek  was  in  roads  or  skid  trails.  The  increases  in 
mean  peak  flow  at  Deer  Creek  also  were  not  significant  (p  =  0.95). 

The  three-day  runoff  volumes  associated  with  high  instantaneous  peak  flows 
were  analyzed.  A  significant  increase  of  1.21  in.  was  detected  on  Needle 
Branch  (Figure  2).  An  increase  of  0.10  in.  for  Deer  Creek  was  not 
significant  (p  =  0.95). 
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Figure    2  — Relation  of  3-day  high-flow   runoff  of  Needle 

Branch  and  of  Deer  Creek  to  that  of  Flynn  Creek.      (Harris, 1977) 
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Harris'  analysis  suggests  that  while  the  average  instantaneous  peak  flow  for 
events  greater  than  50  csm  did  not  increase  significantly,  runoff  volumes 
associated  with  those  flows  did  increase  on  Needle  Branch.  In  effect,  he  is 
suggesting  that  the  duration  of  high  flows  above  some  level  was  increased. 
Harr,  et  al .  (1975)  in  their  analysis  of  peak  flows  at  the  Alsea  Watershed, 
included  smaller  peak  flows  than  those  used  by  Harris.  Their  average  peak 
of  37  csm  was  roughly  one-half  that  of  Harris'  peaks.  What  is  significant, 
however,  is  that  by  including  these  smaller  peaks,  the  average  peak  flow  on 
Needle  Branch  increased  significantly  (p  =  0.99)  by  16.  csm  in  the  fall  and 

10  csm  in  the  winter  after  the  watershed  was  82%  clearcut.  At  Deer  Creek 
(sub-watershed  3)  increases  of  3  csm  in  the  fall  and  12  csm  in  the  winter 
were  both  highly  significant.  Deer  Creek  3  was  65%  clearcut  and  had  12%  of 
the  watershed  in  roads,  landings  and  tractor  skid  trails.  At  Deer  Creek  4, 
fall -period  peak  flows  increased  27  csm.  Fall  peaks  increased  16  csm  and 
winter  period  peaks  increased  10  csm  following  90%  clearcutting  at  Needle 
Branch.  These  increases  were  highly  significant  (p  =  0.99). 

Again,  the  implication  of  Harr  et  al .  is  that  while  increases  in  larger 
peaks  were  difficult  to  detect  (and  apparently  related  to  road  building  and 
soil  compaction),  increases  in  lower  instantaneous  flows  were  significant 
(p  =  0.99).  The  question  becomes  how  much,  if  any,  of  the  increases  in  in- 
stantaneous flows  occurred  as  increases  in  the  duration  of  channel -forming 
flows  (flows  over  50  csm  bankfull  discharge)?  The  next  section  provides 
some  additional  analysis  of  changes  in  the  duration  of  flows  greater  than 
bankfull  in  the  Alsea  Watersheds  following  logging. 

Analysis  of  High-Flow  Duration 

Bankfull  discharge  was  estimated  to  be  50  csm  for  the  Alsea  watersheds. 
The  daily  flows  associated  with  instantaneous  peak  flows  of  50  csm  or 
higher  were  determined  by  regression  analysis  for  each  of  the  three  Alsea 
watersheds.  These  flows  are  shown  below: 

Daily  Flow  Asso- 
ciated with  a  Bank- 
Stream  Bankfull  Flow,  cfs         full  Event,  cfs 

Flynn  Creek  39  30  ( R2  =  0.95) 

Needle  Branch  13  9  (R2  =  0.81) 

Deer  Creek  58  43  (R2  =  0.95) 

Days  when  flow  exceeded  the  threshold  flow  were  then  tallied  for  each 
watershed  and  each  year.  These  results  appear  in  Table  1.  The  numbers 
represent  the  duration  (days)  of  time  that  flow  was  likely  to  be  at  or  above 
bankfull. 

To  test  the  effect  of  timber  harvesting  on  this  variable,  we  did  covariance 
analyses  of  the  high-flow  durations  of  both  Needle  Branch  andJDeer  Creek, 
using  Flynn  Creek  as  the  control  watershed. 
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Table  1.  Duration  of  Flow  at  or  Above  Bankfull  for  Three  AT  sea  Watersheds. 
Water  Year      Flynn  Creek      Needle  Branch      Deer  Creek 


(Days) 

T959 

2 

4 

4 

1960 

2 

2 

2 

1961 

7 

7 

9 

1962 

4 

4 

4 

1963 

2 

3 

3 

1964 

4 

7 

5 

1965 

11 

12 

12 

1966 

8 

9 

9 

1967 

4 

11 

3 

1968 

1 

13 

1 

1969 

7 

15 

6 

1970 

5 

11 

6 

1971 

5 

12 

10 

1972 

8 

19 

9 

1973 

3 

8 

3 

re-treatment  mean 

(1959-65) 

5 

6 

6 

ost-treatment  mean 
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(1967-73) 
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13 
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Figure  3.     Regression  of  Needle  Branch  on  Flynn  Creek  pre-harvest  annual   high 
flow  duration.     Post-harvest  data  points  are  also  plotted  for 
comparison.     The  pre-harvest  relationship  is  (Needle  Branch  =  0.96  X 
[Flynn  Creek]  +  1.20).     R2  =  0.89.  ~  ~ 
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Figure  4.     Regression  of  Deer  Creek  on  Flynn  Creek  pre-harvest  annual 
high  flow  duration.     Post-harvest  data  points  are  also 
plotted  for  comparison.     The  pre-harvest  relationship  is 
(Deer  Creek  =  1 .04  X  [Flynn  Creek]  +  0.80).     R2  =  0.95. 


21 


Clearcutting  of  Needle  Branch  resulted  in  a  123%  average  increase  in  the 
duration  of  high  flows  (p  =  .99).     All   post-treatment  years  showed  increases 
(Fig.  3).     The  analysis  of  Deer  Creek,  which  was  partial-cut,  showed  no 
difference  in  high  flow  duration  as  a  result  of  the  treatment  (p  =  .95). 
One  year,  1971,  had  a  significant  post-treatment  response  (Fig.  4). 

Peak  Flow  Prediction:     Regression  Approach 

Sloan,  et  al .   (1983)  report  on  investigations  by  Corbett  (1979)  on  the 
quickflow  response  of  a  small    forested  watershed  in  Pennsylvania.     The 
antecedent  flow  rate  was  found  to  have  the  greatest  correlation  with 
quickflow  of  all   variables  investigated.     Antecedent  flow  was  believed  to  be 
the  best  overall   indicator  of  the  extent  of  saturated  source  area  within  the 
watershed.     The  purpose  of  this  section  is  to  develop  a  two  parameter  linear 
model   for  peak  flow  prediction  for  the  Alsea  Watersheds  of  the  form 

Qp  =  b0  +  b-\  Qant  +  b2  P24  (6) 

where 

Qp  =  peak  flow  greater  than  50  csm 

Qant  =  mean  daily  streamflow  one  day  prior  to  the  peak  flow 
P24  =  24  hour  precipitation  total   for  the  day  of  peak  flow 
bo>  bl >  b2  =  regression  coefficients 

The  objective  was  to  see  if  (1)  this  model   accurately  predicted  peak  flows 
in  the  Alsea  Watersheds,  and  (2)  any  changes  in  regression  coefficient  from 
the  pre-  to  the  post-logging  period  could  be  attributed  to  logging.     Runoff 
data  is  from  the  USGS  Water  Supply  Papers  for  Oregon.     Precipitation  data  is 
from  the  National   Weather  Service  24-hour  rain  gage  at  Tidewater,  Oregon. 
Ideally,  one-hour  precipitation  data  collected  on  the  watersheds  would  have 
been  available.     Also,  the  Tidewater  gage  was  read  daily  at  5:00  p.m.,   so 
the  24-hour  amounts  don't  correspond  precisely  with  the  24-hour  periods  for 
the  runoff  data. 

Results  of  the  stepwise  multiple  linear  regression  analysis  are  provided  in 
Table  2.     Means  of  variables  are  provided  in  Table  3  and  correlation  coef- 
ficients are  provided  in  Table  4.     Analyses  were  conducted  both  with  all 
flows  greater  than  50  csm  and  with  only  those  events  which  caused  peaks 
greater  than  50  csm  on  Flynn  Creek  (the  "control"  watershed).     The  ratio  of 
the  number  of  50  csm  events  on  Deer  Creek  and  Needle  Branch  to  those  on 
Flynn  Creek  increased  from  1.21:1   and  1.57:1   to  1.38:1   and  2.19:1    for  Deer 
Creek  and  Needle  Branch  respectively  in  the  post-logging  period. 

In  general,  the  regression  model    proved  to  be  a  good  predictor  of  instan- 
taneous peak  flow  based  upon  the  fairly  high  (.6  to  .8)  r*  values. 
However,  the  entire  correlation  structure  within  the  model   changed  in  the 
post-logging  period  -  even  for  the  control  watershed,  so  the  model   lacked 
sensitivity  to  changes  caused  by  management. 
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Stream   Model 


Table  2a 
Predictive  Equations  for  Flows  Over  50  csm  in  Flynn  Creek 


Pre  Logging 
Predictive  Equation 


Post  Logging 
Predictive  Equation 


Flynn  1  variable  Q  =  21.3  +  16.0  P^  .36 

2  variable  Qp  -  -9.4  +  .8  Qant  +  17.6  ?u         .62 

Deer   1  variable  Q  *  22.6  +  17.3  ?u  .53 

2  variable  Q  =  12.3  +  .6  Qant  +  14.6  P^         .77 

Needle  1  variable  Q  ■  19.7  +  22.2  P24  .58 

2  variable  Qant  not  significant 


Qp  ■  -15.1  +  29.9  P24  .67 

Qp  =  -71.0  +  1.4  Qant  +  37.6  P^  .81 

Qp  =  .1  +  26.6  P24  .65 

-  -62.0  +  36.2  P24  +  1.3  Qant  .80 

Qp  =  .6  +  35.9  P24  .69 
Qant  not  significant 


Table  2b 
Predictive  Equations  for  All  Flows  Over  50  csm 


Pre  Logging  Post  Logging 

Predictive  Equation  Predictive  Equation 

Stream   Model  r  -2 


r 


Flynn     1  variable  Qp  -  26.7  +  1.5  Qant  .61  Q     =  -15.1   +  29.9  P^  .67 

2  variable  Qp  =  -34.5  +  1.5  Qant  +  20.0  ?u  .80  Q     =  -71.0  +  1.4  Qant  +  37.6  ?u  .81 

Deer       1  variable  Qp  =  42.4  +  1.2  Qant  .78  Q     =  -2.6  +  25.9  ?u  .59 

2  variable  Qp  =  15.9  +  1.1   Qant  +  9.9  P£4  .85  Q     =  -43.6  +  1.0  Qant  +  32.1   P^  .75 

Needle  1  variable  Q     =  57.0  +  1.1   Qant  .65  Q     =     10.5+30.9  Pp.  .58 

2  variable  Q     =  12.0  +  .9  Qant  +  16.8  P24  .79  Qant  not  significant 


Table  3a 

Means  and  Standard  Deviations  (S.D.)  for  Flows 

Over  50  csm  in  Flynn  Creek 


Stream 


Pre  Logging 

Number 
Qp      Qant    P24   of  Events 


13 


Flynn 

mean 

70.6 

31.0 

3.1 

s.d. 

20.1 

12.3 

.8 

Ui 

Deer 

mean 

75.8 

28.8 

3.1 

s.d. 

18.0 

14.2 

.8 

Needle 

mean 

88.1 

31.3 

3.1 

s.d. 

22.3 

14.7 

.8 

13 


13 


Post 

Logg 

ing 

QP 

Qant 

P24 

Number 
of  Events 

75.1 

23.0 

3.0 

16 

36.9 

11.1 

1.0 

80.3 

25.3 

3.0 

16 

33.1 

12.1 

1.0 

109.0 

35.9 

3.0 

16 

43.7 

16.9 

1.0 

Table  3b 

Means  and  Standard  Deviations  (S.D.) 

for  All   Flows  Over  50  csm 


Pre 

Logging 

Post 

Logging 

Stream 

QP 

Qant 

P24 

Numbi 
of  Evei 

sr 

its 

QP 

Qant 

P24 

Number 
of  Events 

Flynn 

mean 

78.1 

34.6 

3.1 

14 

75.1 

23.0 

3.0 

16 

s.d. 

34.2 

18.0 

0.7 

36.9 

11.1 

1.0 

Deer 

mean 

77.8 

29.9 

3.0 

17 

74.1 

21.4 

3.0 

22 

s.d. 

30.0 

22.4 

0.8 

29.8 

12.5 

0.9 

Needle 

mean 

88.5 

29.8 

3.0 

18 

93.5 

27.2 

2.7 

35 

s.d. 

32.4 

24.6 

0.8 

35.3 

18.0 

0.9 

Table  4a 

Correlations  for  Flows  Over  50  csm 

in  Flynn  Creek 


Pre  Logging 
Correlations 

Stream 

Qp  x  Qant 

Qp  x  p24 

Qant  x 

Flynn 

.43 

.60 

-.12 

Deer 

.64 

.73 

.22 

Needle 

.39 

.76 

.17 

Post  Logging 
Correlations 


Qp  x  Qant      Qp  x  P?4 


Qant  x  P 


24 


07 

.82 

-.49 

19 

.81 

-.61 

26 

.83 

-.23 

Table  4b 
Correlations  for  All  Flows  Over  50  csm 


Pre  Logging 
Correlations 


Post  Logging 
Correlations 


Stream 


Q  x  Qant    Qp  x  P24    Qant  x  P 


24 


Qp  x  Qant 


Qp  x  P24    Qant  x  P 


24 


Flynn 

Deer 

Needle 


.78 
.89 
.80 


.46 
.51 
.63 


.04 
.31 
.33 


-.07 
.03 
.04 


.82 
.77 
.76 


■.49 
.43 
.11 


00 


In  the  pre-logging  period,  for  those  storms  which  caused  peaks  greater  than 
50  csm  on  Flynn  Creek,  antecedent  flow  was  positively  correlated  to  peak 
flow.     Twenty- four  hour  precipitation,  already  the  strongest  predictor 
became  even  more  strongly  positively  correlated  with  post-logging  peaks 
The  ratio  of  Qant/Qpeak  decreased  slightly  from  roughly  0.39:1   to  0.32:1   in 
the  post-logging  period  for  all   three  watersheds. 

When  the  peak-flow  regression  analysis  was  originally  tested,  we  hypothe- 
sized that  a  possible  effect  of  logging  might  be  to  decrease  the  ratio  of 

Qant/Qp  in  the  post-logging  period,  resulting  in  a  higher  peak  flow  for  a 
given  precipitation  event.     This  result  could  also  be  reflected  in  an 
increase  in  the  regression  coefficient  for  P24  which  would  increase  QD 
for  a  given  P24.     While  the  coefficient  for  P24  increased  in  the  post- 
logging  period  it  did  so  on  Flynn  Creek  too.     Therefore,  our  analysis  did 
not  provide  evidence  in  support  of  our  initial   hypotheses.     This  type  of 
analysis  might  be  more  fruitful   given  shorter  rainfall   intensity  data  (e.  g. 
6  hr.),  an  on-site  precipitation  gage  and  a  lower  peak-flow  threshold. 
Hewlett  and  Fortson  (1984)  found,  however,  that  inclusion  of  maximum  hourly 
rainfall   intensity  did  not  greatly  improve  peak  flow  predictions  in  15 
drainage  basins  when  gross  storm  rainfall   and  pre-storm  flow  rate  were 
included  as  independent  variables. 
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RECOMMENDATIONS 

The  following  are  generalized  recommendations,  and  should  be  applied  only  as 
one  aid  to  professional  judgement  and  after  an  understanding  of  their 
background  and  derivations. 

Annual   Flows 

Annual   flows  are  increased  following  logging,  primarily  because  of  reduced 
evapotranspiration.     The  increase  is  dependent  upon  vegetation  type,  annual 
precipitation,  years  since  logging,  and  other  factors.     The  increase  can  be 
estimated  by  equation  (3)  which  was  developed  at  the  H.  J.  Andrews  Forest. 
A  more  site-specific  procedure  recommended     for  estimating  seasonal   or 
annual   runoff  changes  due  to  logging  is  provided  in  WRENS. 

Design  Peak  Flows 

In  low  elevation,  rain  dominated  zones  in  the  Oregon  Coast  Range,  logging 
probably  does  not  have  a  major  effect  on  the  20  to  100-year  return  period 
flows  generally  used  in  culvert  and  bridge  designs.     To  estimate  large 
design  flows  a  Log-Pearson  III  flood  frequency  analysis  with  regionalized 
skew  coefficients  is  recommended  on  gaged  basins.     Any  one  of  several  em- 
pirical relationships,  such  as  equations  (la-Id),  can  be  used  to  estimate 
design  flows  on  ungaged  basins.     We  tried  to  validate  the  utility  of  the  SCS 
Curve  Number  procedure  for  predicting  design  flows,  and  the  effects  (on 
Curve  Number)  of  timber  harvest.     However,  we  could  not  identify  a  suitable 
data  base  for  the  analysis.     Furthermore,  24-hour  design  storms  are  inade- 
quate descriptors  of  major  runoff-producing  events. 

Design  flows  may  be  affected  somewhat  when  over  6-12%  of  a  watershed  is  in 
roads  or  severely  compacted  skid  trails,  landings,  etc.     In  this  case,  a 
small   safety  factor  (15-20%)  may  need  to  be  applied  to  design  flows,  and/or 
extra  considerations  may  need  to  be  given  to  location,  design,  drainage  and 
rehabilitation  measures  for  roads  and  compacted  areas. 

No  consideration  is  given  here  to  peak  flows  affected  by  snow  accumulation 
and  melt. 

Bankfull   Flows 

Clearcut  logging  may  increase  the  duration  of  flows  over  bankfull  by  roughly 
120  percent.     These  increases  will   decay  following  revegetation  and  also 
will   be  proportionately  smaller  as  smaller  percentages  of  the  watershed  are 
logged.     The  effect  of  these  larger  (channel -forming)  flows  on  channel   sta- 
bility and  channel  morphology  will   be  highly  dependent  upon  the  corres- 
ponding increase  in  sediment  delivery  -  particularly  in  the  bed-material 
sizes.     While  channel   stability  may  be  reduced  following  logging  by  the 
increased  duration  of  greater  than  bankfull    flows,  longer-lasting  impacts  on 
channel  morphology  can  be  expected  when  increased  erosion  (particularly  mass 
erosion)  results  in  increased  delivery  of  bed  load  sediment  to- channels. 
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Daily  or  Peak  Discharge  Models 

No  daily  or  peak  discharge  models  were  identified  which  have  been  adequately 
validated  for  use  in  the  Oregon  Coast  Range.     The  Overton-White  Models 
(Overton  and  White,  1978)  and  the  simple  daily  rainfall  model   developed  by 
Sloan,  et  al .  (1983)  deserve  additional  testing,  and  validation  in  the  Coast 
Range.     Further  investigation  of  the  simple  multiple  linear  regression 
approach  as  presented  above  (eqn.  6)  for  predicting  storm-period  runoff 
volume  and  instantaneous  peak  flow  is  also  recommended  -  with  special 
emphasis  on  relating  model   parameters  to  watershed  variables. 


31 


LITERATURE  CITED 

Abdul,  A.  S.  and  R.  W.  Gillham.  1984.  Laboratory  studies  of  the  effects 
of  the  capillary  fringe  on  streamflow  generation.  Water  Resources  Res 
20(6):  691 -698. 

Andrews,  E.D.  1983.  Entrainment  of  gravel  from  naturally  sorted  riverbed 
material.  Geo!.  Soc.  of  Amer.  Bui.  94:1225-1231. 

Beschta,  R.L.  1984.  Long-term  changes  in  channel  widths  of  the  Kowai 
River,  New  Zealand.  Jour.  Hydro!.  (NZ)(In  Press.). 

Brown,  G.W.  1980.  Forestry  and  Water  Quality.  O.S.U.  Bookstores, 
Con/all  is,  Oregon.  124pp. 

Campbell,  A.J.  and  R.C.  Sidle.  1984.  Prediction  of  peak  flows  on  small 
watersheds  in  Oregon  for  use  in  culvert  design.  Water  Resources  Bui. 
20(1):9-14. 

Christner,  J.  and  R.D.  Harr,  1982.  Peak  streamflows  from  the  transcient 
snow  zone,  Western  Cascades,  Oregon,  IN:  Proceedings,  50th  Western  Snow 
Conference,  Reno,  Nevada:  27-38. 

Dunne,  T.  1978.  Field  studies  of  hillslope  flow  processes.  IN:  Kirkby, 
(ed. ),  Hillslope  Hydrology.  J.  Wiley  and  Sons,  New  York,  p. 227-293. 

Dunne,  T.  and  L.B.  Leopold.  1978.  Water  in  Environmental  Planning.  W.  H. 
Freeman  and  Co.,  San  Francisco,  81 8p. 

Emmett,  W.W.  1976.  Bedload  transport  in  two  large,  gravel -bed  rivers, 
Idaho  and  Washington.  Proc.  of  the  Third  Inter-Agency  Sedimentation 
Conference,  Mar.  22-25,  Denver,  CO, :4-101 -4-114. 

Galbraith,  A.F.  1973.  A  water  yield  and  channel  stability  analysis 

procedure.  USDI  Forest  Service,  Kootenai  National  Forest,  Libby,  MT.  38p. 

Goldstein,  R.A. ,  J.B.  Mankin,  and  R.J.  Luxmore.  1974.  Documentation  of 
PROSPER:  a  model  of  atmosphere-soil -pi  ant-water  flow.  Environ.  Sci. 
Div.  Publ.  No.  579.  Oak  Ridge  Natl.  Lab.,  Tenn.  75  p. 

Green,  R.E.  and  J.  C.  Corey.  1971.  Calculation  of  hydraulic  conduc- 
tivity: A  further  evaluation  of  some  predictive  methods.  Proc.  Soil 
Sci.  Soc.  Amer.  35:3-8. 

Harr,  R.D.  1976a.  Forest  practices  and  streamflow  in  western  Oregon. 
USDA,  Forest  Service  General  Technical  Report  PNW-49,  Portland,  Oregon, 
18  p. 


.  1976b.  Hydrology  of  small  forest  streams  in  western  Oregon.  USDA, 

Forest  Service  General  Technical  Report,  PNW-55,  Portland, -Oregon,  15  p. 


32 


.  1980.  Streamf"Iow  after  patch  logging  in  small  drainages  within  the 

Bull  Run  Municipal  Watershed,  Oregon.  USDA  Forest  Service  Research  Paper 
PNW-268,  16  pp. 

.  1982.  Fog  drip  in  the  Bull  Run  Municipal  Watershed,  Oregon.  Water 

Resources  Bui.  18( 5): 785-789. 

.  1983(a).  Effects  of  timber  harvest  on  streamflow  in  the  rain- 
dominated  portion  of  the  Pacific  Northwest.  IN:  Procedings,  Soil,  Water 
and  Fisheries  Technology  Transfer  Symposium,  USDA  Forest  Service,  Pacific 
Northwest  Region,  Portland,  Oregon.  45  pp. 

.  1983(b)  Potential  for  augmenting  water  yield  through  forest 

practices  in  wesern  Washington  and  western  Oregon.  Water  Resources  Bui. 
19(3):  383-393. 

Harr,  R.D.,  W.C.  Harper,  J.  T.  Krygier,  and  F.  S.  Hsieh.  1975.  Changes 
in  storm  hydrographs  after  road  building  and  clear-cutting  in  the  Oregon 
Coast  Range.  Water  Resources  Res.  II( 3): 436-444. 

Harr,  R.D.,  and  R.M.  McCorison.  1979.  Initial  effects  of  clearcut  log- 
ging on  size  and  timing  of  peak  flows  in  a  small  watershed  in  western 
Oregon.  Water  Resour.  Res.  15,  p.  90-94. 

Harr,  R.D.,  A.  Levno,  and  R.  Mersereau,  1982.  Streamflow  changes  after 
logging  130-year-old  Douglas  fir  in  two  small  watersheds.  Water 
Resources  Res.  18( 3): 637-644. 

Harr,  R.D.  and  S.N.  Berris,  1983.  Snow  accumulation  and  subsequent  melt 
during  rainfall  in  forested  and  clearcut  plots  in  western  Oregon.  IN. 
Proceedings,  51st  Western  Snow  Conference,  Vancouver,  Washington:  33P45. 

Harris,  D.D.  1977.  Hydro! ogic  changes  after  logging  in  two  small  Oregon 
coastal  watersheds.  USDI  Geo!.  Surv.  Water  Supply  Pap.  2037.  31  pp. 

Harris,  D.D.,  L.L.  Hubbard,  and  L.  E.  Hubbard,  1979.  Magnitude  and  fre- 
quency of  floods  in  western  Oregon.  USGS  Open  File  Report  79-553, 
Portland,  Oregon.  35  pp. 

Hewlett,  J.D.  and  J.C.  Fortson,  1984.  Additional  tests  on  the  effect  of 
rainfall  intensity  on  storm  flow  and  peak  flow  from  wild-land  basins. 
Water  Resources  Res.  20(7),  985-989. 

Hewlett,  J.D.  and  A.R.  Hibbert.  1967.  Factors  affecting  the  response  of 
small  watersheds  to  precipitation  in  humid  areas,  IN:  W.  E.  Sopper  and  H. 
W.  Lull  (Leds.).  Int.  Symp.  for  Hydro!.,  Pergamon  Press,  New  York,  New 
York,  p.  275-290. 

Jackson,  W.L.  and  R.L.  Beschta,  1982.  A  model  of  two-phase  -bedl-oad 

transport  in  an  Oregon  Coast  Range  stream.  Earth  Surface  Processes  and 
Landforms,  Vol.  7,  p.  517-527. 


33 


Jackson,  W.L.   and  R.L.   Beschta,   1984.     Influences  of  increased  sand 

delivery  on  the  morphology  of  sand  and  gravel   channels.     Water  Resourdes 
Bui.   (IN  PRESS). 

Johnson,  M.G.   and  R.L.   Beschta.     1980.     Logging,   infiltration  capacity, 
and  surface  erodibility  in  western  Oregon.     Jnl .   of  Forestry,  June, 
1980.     334-337. 

Leaf,  C.F.  and  G.E.   Brink.     1973.     Hydrologic  simulation  model   of  Colorado 
subalpine  forest.     USDA,   forest  Service  Res.   Paper  RM-107,  23  pp. 


.     1975.     Land  use  simulation  model   of  the  subalpine  coniferous  forest 

zone.     USDA  Forest  Service  Res.   Paper  RM-135,  42  pp. 

Leavesley,  G.H.,  R.W.   Lichty,  B.M.Troutman  and  L.G.   Saindon.     1983.     Pre- 
cipitation-Runoff modeling  system:   users  manual.     USDI  Geological   Survey, 
Water  Resources  Investigations  Report  83-4238,  Denver,  Colorado.     207  pp. 

Lystrom,  D.J.     1970.     Evaluation  of  streamf low-data  program  in  Oregon. 
USDI  Geo!.  Surv.   Open  File  Report.     Water  Resources  Division.     Portland, 
OR     .28  pp. 

Leopold,  L.B.   and  T.  Maddock  Jr.     1953.     The  hydraulic  geometry  of  stream 
channels  and  some  physiographic  implications.     USDI  Geo!.   Surv.   Prof. 
Paper  252,   57  pp. 

Lyons,  J.K.   and  R.L.   Beschta.     1983.     Land  use,   floods  and  channel   changes: 
Upper  Middle  Fork  Willamette  River,  Oregon  ( L1936-1980). 

Moring,  J.R.   and  R.L.   Lantz.     1975.     The  Alsea  Watershed  Study:     Part  I  - 
Biological  Studies.     Fisheries  Research  Report  Mo.   9,  Oregon  Dept.  Fish 
and  Wildlife,  Corvallis,  Oregon.   66  pp. 

Osborn,  J.F.     1980.     Estimating  streamflow  characteristics  at  spawning  sites 
in  Oregon,  CERL-051 ,  U.S.   Environmental   Protection  Agency,  Corvallis, 
Oregon,  96  pp. 

Overton,  W.S.   and  C.   White.     1978.     Evolution  of  a  hydrology  model   -  an 
exercise  in  modelling  strategy.     Int.  J.   General   Systems,  Vol.  4:89-104. 

Reid,  L.M.   and  T.   Dunne.     1984.     Sediment  production  from  forest  road 
surfaces.     Water  Resources  Res.   20(11),   1753-1761. 

Rothacher,  J.     1963.     Net  precipitation  under  a  Douglas  fir  forest.     For. 
Sci.    9(4): 423-429. 


1971.     Regimes  of  streamflow  and  their  modification  by  logging. 
IN:  J.  Morris  (ed.),  Forest  land  uses  and  stream  environment,  symposium 
proceedings,   1970,   Oregon  State  University,   Corvallis,   OR,   pp  40-54. 


34 


Schumm,  S.A.  1971.  Fluvial  Geomorphology:  channel  adjustment  and  river 
metamorphosis,  IN:  H.W.  Shen  (ed. ),  River  Mechanics,  Yol .  I,  Ft.  Collins 
Coloradfo.  5.1-5.22. 

Silvey,  L.  and  D.  Rosgen.  1981.  HYSED:  A  water  resource  analysis  model 
for  forest  land  use  planning.  USDA  Forest  Service,  Rocky  Mountain 
Region.  Lakewood,  Colorado.  62  p. 

Simons,  D.B.,  R.M.  Li  and  M.A.  Stevens,  1975.  Development  of  models  for 
predicting  water  and  sediment  routing  and  yield  from  storms  on  small 
watersheds.  Engineering  Research  Center,  Colorado  State  University  Ft 
Collins,  Colorado,  130  pp. 

Sloan,  P.G.,  I.D.  Moore,  G.B.  Coltharp,  J.D.  Eigel.  1983.  Modeling 

Surface  and  Subsurface  Stormflow  on  Steeply-Sloping  Forested  Watersheds. 
Research  Report  142,  Water  Resources  Res.  Inst.  Univ.  of  Kentucky, 
Lexington,  Kentucky.  157  pp. 

Swanson,  F.J.  and  G.W.  Lienkaemper,  1978.  Physical  consequences  of  large 
organic  debris  in  pacific  northwest  streams.  USDA  Forest  Service  General 
Technical  Report  PNW-69,  12  pp. 

Troendle,  C.A.  and  J.D.  Hewlett.  1979.  A  variable  source  area  hydrograph 
simulator  (VSAS)  for  small  forested  watersheds.  Unpub.  paper,  Univ.  of 
Georgia,  Athens,  Georgia,  38  pp. 

Troendle,  C.A.  and  C.F.  Leaf.  1980.  Hydrology,  IN:  An  approach  to  water 
resources  evaluation  of  non-point  silviculturaT~sources.  USDA,  Forest 
Service.  EPA-600/8-80-012,  contact  EPA-IAG-DG-0660,  Chapter  3. 


1 


35 


I 


APPENDIX  I 

Prosper  Model  j. 

Overton  -  White  Models  I- 

Sloan,  et  al .  Daily  Runoff  Model  I- 
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Figure  111. B.2.— Schematic  of  PROSPER  (from  Goldstein  and  others  1974). 
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i     FLEXFORM  ABSTRACT     Hydrology  Model  II 


Resolution,  one  day; 
units,  m3/ha 
(Figure  2-b) 


Date:  4/24/73 
Revised:  7/11/73 


Model  11-a 
State  Variables 


X- 
*3 

X4 
.X, 


Canopy  storage 
Snow  storage 
Soil  water 
Ground  water 
Streamflow  (dailv 


Input  Variables 


-A.  OoIT| 


h 


Precipitation  as  rain 
Precipitation  as  snow 
Thrufall.     ll-b9)  =  Interception     para- 
meter 

Potential  Evapotranspiration  (PET) 
Average  temperature  (°F) 


A.    Canopy  Storage  and  Evapotranspiration  (ET) 

gj    =(bi0--X])(l-exp(-inr,))  Input  to  x„ 
where  b10  =  maximum  storage  capacity, 

2>u  =  (l-M/bvo. 
g2   -mK{[«i  +  l]-£ai+l]*,'*,,,0}    Adjusted 
PET,  see  Figure  6. 


eg 

en 


'\Z 


FIGURE  6     Adjusted  atmospheric  demand  function  for  Model 

n. 

g3  -niin{jf,+gtl  g2(l  -exp(-fc13U-i  +gi )))} 
Canopy  evaporation,  where  b13 
=  evaDoration  rate. 
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g4  =  Transpiration,  see  Figure  7.  fc,=  wilting 
point.  bz  =  transpiration  resistance  point. 


FIGURE  7     Actual  transpiration  function  for  Model  II. 
g5   =  (1  —bq)il  —  gj  Drip  to  forest  floor. 

B.  Snow  Melt 

g6    =max  ;0.(s3-32)(254RAD  +  0.014(_-4 

+  g5 ))}  Potential  snow  melt,  where  RAD 
is  a  monthly  radiation  index. 

g7   =  min;g6.  r2+x2]  Actual  snow  melt. 

C.  Infiltration 

gi  2  =  =4.  +  Is  -r  g?  Potential  infiltration. 


S13 


=  min  fc->  —  x 


Si  2.1 


Actual    infiltration. 


where  fr7  =  maximum  x3. 
gij.  =  max{gl2  —  gi3.0)  Surface  runoff. 

D.  Soil  Water  Flow 

g8  =  max{<i>1(x3  —  b5),  0}  Total  flow,  where  b5 
^retention  capacity.  b3  =  instantaneous 
flow  rate.  <p1  =  1  -exp(  —  r>3)  =  discrete 
flow  rate. 

g9  =  min{68  — x4,  g8}  Percolation,  where  bs 
=  maximum  x4.  • 


?,^  =   Po 


gQ  Lateral  flow. 


E.    Ground  Water  Lateral  Flow 

gn  =  max{<p2(x4  — £>6),  0}  Total  flow,  where  b6 
=  retention  capacity,  fe4  =  instantaneous 
flow  rate,  <i2  =  l  —  exp(  — b4  )  =  discrete 
flow  rate. 


F.    Stream  Flow 


'Su+Sia.  Stream  flow. 


Model  11-b 


Changed  value  for  bx  of  function  g4  to  see  its  effe 
on  transpiration. 

Model  Il-c 

Readjusted  parameter  values  and  initial  conditions 
and  used  adjusted  driving  data  for  Watershed 
Aiso  used  this  form  with  Watershed  10  driving  dai 
This  form  may  be  direct!}   compared  with  iai 
model*,  while  Il-a  and  l!-b  ma^  not, 


FLEXFORM  ABSTRACT   Hydrology  Model  I 
Resolution,  one  day:  units,  m3  Date:  5/4/7 

(Figure  2-c.  d) 


Model  Ill-a 

State  Variables 

As  in  Model  II,  plus 

x22     Lagged  ground  water 

Input  Variables 

As  in  Model  II,  except 

zA  =  b1z1     Thrufall,  where   (1  —  c2)  =  Intercept!  < 

parameter 

A.    Canopy  Storage  and  ET 

Same  as  Model  II,  except  for  units  adjustment 


and  notational  changes. 


B.  Snow  Melt 

Same  as  Model  II,  except  for  units  adjustmen 
and  notational  changes. 

C.  Infiltration 

g7   =z4  +  g4  +  g6  Potential  infiltration,  where 
g4  =  drip,  g6  =  actual  snow  melt. 

g14  =  max{0,  (x3-fcu)  +  g?  +  (g11-g8)-gi: 
—  g13}  Surface  runoff,  where  bu 
=  maximum  x3,  g13=  transpiration;  g 
gn,gi2  see  below. 

g7-g14  Actual  infiltration. 
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D.  Soil  Water  Flow 


g6  =max{0,  </>i((x3-b9)  +  blsg-)}  Percol- 
ation, where  pretention  capacity,  b-> 
=  instantaneous  flow  rate,  <b  =1—  exp 
(— b7)  =  discrete  flow  rate,  b,s 
=  "residence"  parameter'  which  is  the 
proportion  of  g7  which  contributes  to  flow 
during  the  same  day. 


_  jgii 


>     gll  ^OH| 

Ui +*i*» (gn-g,),  gn>g8J 

Lateral    flow,    where    b.s  =  "residence" 
parameter  for  (gll~gs).^l  as  in  gg. 

No  water  moves  laterally  in  soil  until  .v,  is 
filled. 

E.    Ground  Water  Lateral  Flow 

g9  =b17(x22  +  b16g8)  Current  lag  effect,  where 
ij 6  =  "residence"  parameter  for  g.    fc,, 

t  r  DO'       "1    / 

=  lag  factor. 

f22.2i  =  g9—  x22-  Implies  that  current  g9  be- 
comes x22  for  next  time  step. 

g10  =  max{0,  <fi2((x*-b10)  +  bi6ga-g9)}  Lat- 
eral flow,  where  b10=retention  capacity. 
bg=  instantaneous  flow  rate,  <£,  =  1 1' 
exp  (  —  b8)~ discrete  flow  rate. 

g, i=  max {0AxA~b12)  +  ga-glo\  Back  flow  of 
excess  over  storage  capacity,  where  bl2 
=  maximum  x4.  Necessary  device  in  stra- 
tified model  since  uphill  strata  ground 
water  lateral  flow  can  exceed  storage 
capacity  of  downhill  strata. 

F.    Stream  Flow 

gi5  =  gio  +  gu  +  g14  Stream  flow. 

Model  111-b 

Incorporates  a  split  ground  water  compartment 
channelized  (xA)  and  non-channelized  (x23),  im- 
posed on  Model  III-a.  Storage  above  retention 
capacity  is  split  between  the  two  with  x23  receiving 
3°%  (b20=0.30).  Storage  below  retention  capacity 
is  assigned  to  x4;  assignment  to  x23  would  not 
change  behavior. 

Percolation  is  split  with  x23  receiving  50%  (b19 
=  0.50).  Excess  over  x23  storage  capacity  is  sent  to 
xA  where  it  contributes  to  flow  (or  back'flow)  from 
there.  Channelized  lateral  flow  rate  (cp2)  is  left 
unchanged.  Non-channelized  lateral  flow  is  con- 


siderably slower,  but  of  the  same  form,  and  flows 
directly  into  the  stream. 

Model  IH-c 

Identical  to  Ill-b  except  that  a  portion,  bZi  =0.03, 
of  calculated  percolation  passed  directly"  to  the 
stream.  b22  =0  in  subsequent  Model  III  forms,  but 
was  used  again  for  Model  V  only. 

Model  III-d  (with  "dry  sponge") 

The  retention  capacity  ol  x3  was  made  a  function  of 
.v,,  so  that  water  entering  a  nearly  dry  compartment 
would  pass  through  quickly. 

■fh:i-x3<bi4 

gi9  =  <  -v3,  blA£xzSb9 

x  t  bq,  x3>b9 

Retention  capacity  of  xs,  where  b9 
=  retention  capacity  when  wet,  blA 
=  transpiration  resistance  point. 

g16  =  max{0,cp}((x3-Sl9)  +  blsg1)}  Total  per- 
colation. 

gs   =gi6-gi7  Percolation  into  x4,  where  g]7 
=  percolation  into  x13. 

Model  Ill-e 

Model  III-d  with  -2=0,  so  that  no  snow  occurs 
regardless  of  temperature. 


FLEXFORM  ABSTRACT  Hydrology  Model  VII 
Resolution,  one  day;  1  Date- 12/20/73 

u,?its'm3  Revised:  10/75 

(Figure  2-h) 


Model  VII-a 
State  Variables 


X2 

X3 

X4 

xs 

xl 

xs 
x9 
x29 


Foliage  storage 

Epiphyte  storage 

Snow  pack  storage 

Potential  infiltration 

Channelized  soil  water 

Non-channelized  soil  water 

Channelized  ground  water 

Non-channelized  ground  water 

Stream  flow  (daily) 

Lagged  channelized  ground  water 


L 
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Input  Variables 

Average  temperature  (°F) 
Precipitation  as  rain 
Precipitation  as  snow 
Thrufall 
PET 


g22=max 


?12  +  gl4 


-^13(^17-^15)- 
gl5+g20-g21 


A.  Canopy  Storage  and  ET 

gj  =  (l  —  b2)bi/b3bA  Rate  of  Xj  change,  where 
(]  —£■,)  =  interception  parameter.  b5 
=  proportion  of  interception  taken  by 
foliage,  b-  —  maximum  total  canopy  stor- 
age, bj.  =« proportion  of  total  in  foliage 
storage. 

g:  =  (i  -ti2  )(i  —bi) ,.'fcs(!  -bA)  Rate  of  x2 
change. 

g3   =(fc3i>it-x1)(l-exp(-g1z2))  Addition  to 

Xj. 

g4  =(fc3(l-^)-x2)(l-exp(-g2z2))  Ad- 
dition to  x2. 

gs   m*2~ zs  ~ S3  - g4  Canopy  drip. 

g6  =max{0,  [r6  +  l]-[--6  +  l]--^}  Adjusted 
PET.  where  bh  =  rainfall  needed  to  stop 
evaporation  (graph  similar  to  Figure  6). 

g7  =min{.\-1-rg3.  g6(l-exp(-b7U'i+g3)))] 
Evaporation  from  xu  where  b-  =Xj  evap- 
oration rate. 

g8  =  min'b5.(g6-g-),  x2  +  g4}  Evaporation 
from  x2,  where  b8  —  proportion  of  remain- 
ing ET  demand  satisfied  by  x2. 

g13  =  Transpiration,  see  Figure  7,  but  with 
notational  changes. 

B.  Snow  Melt 

Potential  (g9)  and  actual  (g10)  snow  melt  as  in 
Model  III.  with  some  notational  changes. 

C.  Infiltration 

gi  1  —  -5  +  gs  +  gi  0  Potential  infiltration. 


(l-b21) 

iiT8i:T«»~ SUTS20- 521/  Surface 
runoff,  where  £>13(b17  —  bi5)  =  amount  of 
soil  water  storage  above  retention  cap- 
acity which  is  channelized,  (1  —  &21)gu 
=  infiltration  into  x5;  g12,  g14,  g15,  g20, 
g2i  see  below. 

•  g22  Actual  infiltration. 


D.  and  E.  Soil  Water  Flow  and  Ground  Water 
Lateral  Flow 

Percolation,  lateral  flow,  and  back  flow  lof  execs- 
over  storage  capacity)  are  all  modelled  as  estab- 
lished in  III— a;  percolation  to  channelized  ground 
water  is  lagged  with  regard  to  lateral  flow.  The  order 
of  flow  operation  is  (see  Figure  4): 


512 
gl4 

git 

gib 


S19 

g20 
g21 


Lateral  flow  from  x6  to  x5 ; 

Surface  flow  from  x6  to  x5 ; 

Percolation  from  xs; 

Percolation  from  x5  to  x8  (remainder  to 

Lateral  flow  from  x8  to  x7; 

Amount  of  water  currently  withheld  from 

contributing  to  x7  lateral  flow; 

Lateral  flow  from  x7  to  x9 ; 

Back  flow  from  x7  to  x3 ; 

Lateral  flow  from  x5  to  x9. 


F.    Stream  Flow 

g23  =  gi9  +  g2i  +g22  Stream  flow. 

Model  Vll-b  \ 

E.    Ground  Water  Lateral  Flow 

The  lag  function  was  redefined  so  that  it  was  as  first 
conceptualized.  This  removed  the  multiplication  by 
£>2S,  the  "residence'"  parameter  for  percolation  into 
x7,  during  the  updating  of  x29. 

For  biographies  and  photographs  of  the  authors,  please  see  p. 
137  of  this  issue. 


SLOAN    ,et.al. 
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Figure  4.1     Schematic  Flow   Diagram   of  the   Daily 
Watershed  Model. 
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Table  4.1     Watershed  Model  Function  Descriptions 


FUNCTION 


EQUATION 


FUNCTION    SCHEMATIC        PROCESS 


CMAX  =  CEPMAX*  FCAN 


INTERCEPTION 


TIME    OF    YEAR 


RUNOI  =   PB    X     PREC1P 
INFIL  =  (l-PB)    *     PRECIP 


PB  =  FSTP+  PCe 


PAC  a  (US1N/USMAX) 


VARIABLE 
SOURCE  AREA 
RUNOFF 


USWP  USMAX 

SOIL    ZONE    (USIN) 


RUN02  =  Kl  X   FFU 
PERCO  =  (l-K!)  X  FFU 


FFU  =  FU  * 


[uSMAXJ 


SOIL  ZONE 
DRAINAGE    3 
INTERFLOW 


USWP  USMAX 

SOIL  ZONE     (  USIN) 


AEVAP  =   EVAP  (EVAP  <  PE) 

8  PE   (EVAP>  PE) 
EVAP=  (USIN -USWP) 
ERATE 


/  PE  =  7mm/d 
PE  =  4mm/d 


SOIL    ZONE 
EVAPOTRANSPIRAT 


USWP  USMAX 

SOIL  ZONE  (USIN  ) 


RUN03  =  K2  *  FFS 
GW  =  (I-K2)*   FFS 


FFS  s  FSX  (SSIN) 


KS 


~io 

I 

£    6 

S^ 

w    4 

^jt< 

09 

u.    2 

j/T                                                                               I 

u.     n 

GROUNDWATER 
ZONE  SEEPAGE 
a  GROUNDWATER 
FLOW 


GROUNDWATER  ZONE 
(SSIN) 


-72- 


.  "  -■:■"  :.■■  ".■  :,_..■.:■■.  ':...-. ■r-.:/-'^:  :■■■  ■:■•.:■■'■'-  .". 


i&g&g^ 


■B 


Table  4.2     Model  Parameter   Descriptions  and  Values 


Process/Zone 


Parameter 


Definition 


Parameter   Vaiue 
(Little   Millseat    Watershed) 


interception 


Evapotranspi  ra- 
tion 


CEPMAX 
FCAN 


Variable  Source 

FSTR 

Area   Runoff 

PAC 

PC 

Soil   Zone 

USMAX 

KU 

FU 
Kl 

USWP 
ERATE 


Groundwater  Zone      FS 

KS 
K2 


OTHER  VARIABLES 


CM  AX 
USIN 
SSIN 
PB 


Maximum    interception 
capacity  (mm) 

Canopy  development   func- 
tion:   modifies   CEPMAX    for 
time   of   year   (i.e.   canopy- 
development) 

Fraction   of    watershed  always 
contributing   to  direct    runoff 
(i.e.  area   of   stream) 

Source  area  exponent 

Source  area  coefficient 

Soil  zone   thickness   (mm) 

Soil    water   conductivity  ex- 
ponent  (KU=2b+3,  where   -b 
is   the  slope  of   a   log-log 
plot   of   the  soil   water    re- 
tention curve) 

Soil    water  conductivity   coef- 
ficient 

Fraction  of   Soil  Zone  drain- 
age  becoming   interflow 

Wilting  point  water  content 
(input  as  %  by  volume,  used 
as   mm   of   water  in  program) 

Evapotranspiration   rate  coef- 
ficient 

i 

Groundwater  exponent   (1   for 
linear   groundwater   recession) 

Groundwater   recession  constant 

Fraction   of   groundwater   drain- 
age becoming  baseflow 


Actual   interception  capacity  (mm) 
Actual   soil   water   volume   (mm) 
Actual   groundwater  volume   (mm) 
Fraction  of   water  contributing 
to  direct    runoff 


2.02 
See  Table   4.1 

0.05    (0.05) 


39.295    (40)** 

4.11xl0~6 
(4.1   x   10-6)" 
1087   (1070) 

11.810  (1  1.467) 


1.49x10 
1.0   (1.0) 

124  (130) 
11.44%   (12.14%) 

27.4 


*      Groundwater  Zone   does 

Values  in  parentheses  a 

Values  used  in  BROOK 
Watershed 


not   exist   in  the   Little   Millseat   watershed. 
re   the  initial   parameter  estimates  prior   to  optimization 
mode!   (Federer  and   Lash,   1978)   for   Hubbard   Brook 
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